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Abstract 

In this letter the exphcit form of evolution operator for the four atoms Tavis- 
'^~>! ■ Cummings model is given. 

00 ■ This paper is a sequel to the papers [T] and the purpose is to give an exphcit form to 

^_. . the evolution operator of Tavis-Cummings model (j2j) with four atoms. This model is a 

O ■ very important one in Quantum Optics and has been studied widely, see jHj as general 

^ I textbooks in quantum optics. 

We are studying a quantum computation and therefore want to study the model from 

" p ■ this point of view, namely the quantum computation based on atoms of laser-cooled and 

trapped linearly in a cavity. We must in this model construct a controlled NOT gate or 

H ! other controlled unitary gates to perform a quantum computation, see jl] as a general 

^ I introduction to this subject. 

^1 For that aim we need the explicit form of evolution operator of the models with one, 

two and three atoms (at least). As to the model of one atom or two atoms it is more or 
S^ . less known (see [5J), while it was recently given by us j3] for the three atoms. 

Since we moreover succeeded in finding the explicit form for the case of four atoms we 
report it. 

The Tavis-Cummings model (with n-atoms) that we will treat in this paper can be 
written as follows (we set h = 1 for simplicity). 

An n 
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where uj is the frequency of radiation field, A the energy difference of two level atoms, 
a and a^ are annihilation and creation operators of the field, and g a coupling constant, 

'i 1 ^i 



and L = 2". Here en \ en and aj are given as 



of = I2 ® • • ■ ® I2 ® cr^ (S) I2 ® ■ ■ ■ ® I2 {i- position) G M(L, C) (2) 

where s is +, — and 3 respectively and 

^+=(0 oj' ^- = (1 oj' ^'=(0 -ij' ^2=(o ij- ^2) 

Here let us rewrite the hamiltonian (^ . If we set 

n n 1 " 

s^ = Y.^^ s- = T.^t\ 53 = ^E^f^ (4) 

then (PJ can be written as 

H = uj1l® a^a + AS3 ® 1 + ^ (5+ ® a + 5_ ® a"^) =Ho + V, (5) 

which is very clear. We note that {S^, S^, S3} satisfy the s'u(2)-relation 

[S'3, S+] = S+, [S3, S^] = -S_, [S+, S'_] = 2S3. (6) 

However, the representation p defined by p(cr+) = S+, p{cr_) = S^, ^(0-3/2) = 5*3 is a 
reducible representation of sm(2). 

We would like to solve the Schrodinger equation 

Zj^U = HU= {Ho + V) U, (7) 

where U is a. unitary operator (called the evolution operator). We can solve this equation 
by using the method of constant variation. The result is well-known to be 

under the resonance condition A = u, where we have dropped the constant unitary 
operator for simplicity. Therefore we have only to calculate the term (jHJ explicitly, which 
is however a very hard task ^. In the following we set 

An = S+0a + S_^a'^ (9) 

for simplicity. We can determine e"**^"^" for n = 4 (four atoms case) completely. As to 
the cases of n = 1 ~ 3 see [l] . 



^the situation is very similar to that of the paper quant-ph/0312060 in [7| 



Four Atoms Case In this case A4 in (jH)) is written as 
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If we set T as 
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(10) 
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then it is not difficult to see 
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This means a well-known decomposition of spin 



1111, ,11 

-®-®-®- = (0©l)(g)-(g)- 
2 2 2 2^ ^22 



1 1 











1©2. 



Since we have calculated e "^^^^^ in |ilj we have only to do e '^^^^^ which is however very 
hard. The result is 

exp {-itgB2) = 

( f2{N + 2) hi{N + 2)a^ ko{N + 2)a^ \ 

/i(A^+l) ho{N + l)a^ 

/ii(iV)(at)2 foiN) h^i{N)a^ 

/io(A^-l)(at)2 /„i(Ar-l) 

V A;o(iV-2)(at)4 /i_i(iV - 2)(at)2 /„2(iV-2)/ 

/ -2zFi(A^ + 2)a -2iHo{N + 2)a^ \ 

-2iFi(iV + l)at -iHi{N+l)a -2iifo(iV + l)a^ 

-lHi{N)a^ -i/7_i(A^)a 



+ 



-2iHo{N - l){a^)'^ 







V 







-2«i7o(A^-2)(at) 



t^3 



\H^i{N-l)a^ 









-2iF_i(A^-2)a 



-2iF^i{N -l)a 




'11^ 



where 



/2 = 1 + 4(iV- l){(n+/A+)(cost^^- 1) - (M_/A_)(cost(7yAr- l)}/v^ 

/i = («+ costgJ\+ — U- costgJ\^)/vd, 

fo = l + 2{{v+w+/ X+) {cos tg^ - 1) - (i;_w_/A_)(cost^^ - 1)}/Vd, 



f_i = (m_|_ costgJX_ — u_ costgJX^)/vd, 

f.2 = l + 4(iV + 2){{u+/X^){costg^ - 1) - {u_/X+){costg^ - 1)}/Vd, 
h = 2{{v+/X+){costg^-l) - {v^/X^){costg^ - 1)}/Vd, 
ho = {costgJX^ — costgJX^)/Vd, 

h_i = 2{{w+/ X+) {cos tg ,/x^ - 1) - {w_/X^){costg,/x^ - 1)}/Vd, 
fco = 4{{l/X+){costg^-l) - {l/X^){costg^-l)}/Vd, 
and 

^1 = {(^+/\A+) sint^-yA^- (^-/\A^) smtg^Jx^}/Vd, 
F-i = {{u+/^X+) sintg^/xl - {u_/^xl) smtg^xl} /\/d, 
Hi = 2{{v+/^+) smtg^-{v^/^) smtg^}/Vd, 
Ho = {{l/^) sintg^ - (1/^) smtg^}/Vd, 
i7_i = 2{{w+/^+) smtg^+ - {w^/^-) smtg^}/Vd 
, and d = d{N), X± = X±{N), u± = u±{N), v± = v±{N) and w± = w±{N) defined by 

X±{N) = 10N + 5± 3^d(N), u±{N) = ^(-3 ± \fd(N)), 

v±iN) ^^^i2N-l± yd(AO), w^{N) ^ \l^i2N + 3 ± ^/d{N)), 
d{N) = AN^ + 4A^ + 9. 



This form is very complicated. We note tliat in tlie process of calculation we used 
Mathematica to the fullest. 

A comment is in order. We would like to generalize the results in this paper and P 
to the cases of more than four atoms. However, it is not easy to perform it due to some 
technical reasons. There is a (big ?) gap between the four atoms and the five ones. 



We obtained the explicit form of evolution operator of the Tavis-Cummings model for 
the case of three and four atoms, so many applications to quantum optics or mathematical 
physics will be expected, see for example . In the near future we will apply the result to 
a quantum computation based on atoms of laser-cooled and trapped linearly in a cavity 

We conclude this paper by making a comment. The Tavis-Cummings model is based 
on (only) two energy levels of atoms. However, an atom has in general infinitely many 
energy levels, so it is natural to use this possibility. We are also studying a quantum 
computation based on multi-level systems of atoms (a qudit theory) j7j. Therefore we 
would like to extend the Tavis-Cummings model based on two-levels to a model based 
on multi-levels. This is a very challenging task. 
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